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ABSTRACT 

Radio pulsar PSR J1028-5819 was recently discovered in a high-frequency search (at 3.1 GHz) in the 
error circle of the EGRET source 3EG J1027-5817. The spin-down power of this young pulsar is great 
enough to make it very likely the counterpart for the EGRET source. We report here the discovery of 
7-ray pulsations from PSR J1028-5819 in early observations by the Large Area Telescope (LAT) on 
the Fermi Gamma-Ray Space Telescope. The 7-ray light curve shows two sharp peaks having phase 
separation of 0.460 ±0.004, trailing the very narrow radio pulse by 0.200 ±0.003 in phase, very similar 
to that of other known 7-ray pulsars. The measured 7-ray flux gives an efficiency for the pulsar of 
~ 10 — 20% (for outer magnetosphere beam models). No evidence of a surrounding pulsar wind nebula 
is seen in the current Fermi data but limits on associated emission are weak because the source lies in 
a crowded region with high background emission. However, the improved angular resolution afforded 
by the LAT enables the disentanglement of the previous COS-B and EGRET source detections into 
at least two distinct sources, one of which is now identified as PSR J1028-5819. 
Subject headings: pulsars: general — stars: neutron 
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1. INTRODUCTION 

One of the most intriguing legacies of the Energetic 

Gamma-Ray Experiment Telescope (EGRET) on the 
Compton Gamma- Ray Observatory (CGRO) was the 
group of 150 sources (Hartman ct al. 1999) that could 
not be firmly identified with any known counterparts. 
The majority of these (about 100) are in a population 
clustered toward the Galactic plane, and many were pre- 
sumed to be 7-ray pulsars. Indeed, a number of radio 
pulsars were discovered to lie in or near EGRET uniden- 
tified source error boxes after the end of CGRO obser- 
vations in 2000, both in large surveys such as the Parkes 
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Multibeam Survey (Kramer et al. 2003) or in deep obser- 
vations in some of the EGRET error boxes (e.g. Roberts 
et al. 2002, Halpernet al. 2001). Unfortunately, searches 
for pulsations at the radio periods in the EGRET archival 
data were not feasible, given the small number of 7-ray 
photons in these sources and the difficulty of predicting 
pulsar phase in the presence of rotational instabilities 
(also known as timing noise). 

The situation has improved dramatically with the re- 
cent launch of Fermi on June 11, 2008, which has been 
operating successfully through early calibration and now 
in sky survey mode. The Fermi LAT has a sensitivity 
that is more than an order of magnitude that of EGRET, 
and it is already possible to better perform searches for 7- 
ray pulsations in many of the EGRET sources. This Let- 
ter reports results on one such search in the error circle 
of the EGRET source 3EG J1027-5817, at the location 
of the radio pulsar PSR J1028-5819, discovered (Keith et 
al. 2008) just a few months prior to the launch of Fermi 
as part of a search of three EGRET sources at high fre- 
quency. PSR J1028-5819 is a young pulsar, with period 
P = 91.4 ms, period derivative P = 1.61 x 10^^** ss^^ 
and characteristic age of 9.21 x 10^ yr. The derived spin- 
down power, = 8.43 x lO^^ergs"^ combined with its 
dispersion measure-derived distance of 2.3 kpc makes it a 
plausible counterpart for the EGRET source with flux of 
6.6 ±0.7 X lO-^phcm-^s"^ {E > 100 MeV). The radio 
pulse profile is extremely narrow, consisting of two highly 
linearly polarized components. The full width is only 560 
/is, giving it the smallest duty cycle of any known pul- 
sar and an order of magnitude smaller duty cycle than 
other young pulsars (e.g. Weltevrede & Johnston 2008). 
Fermi LAT has now discovered the 7-ray pulsations from 
PSR J1028-5819, and we can thus confirm that some of 
the photons attributed to 3EG J1027-5817 originate from 
PSR J1028-5819. 

2. OBSERVATIONS 

Fermi was launched into low-Earth orbit and, after 
a six-week commissioning phase, began nominal sky- 
survey observations on August 11, 2008. The LAT, the 
main instrument on Fermi, is a pair-production telescope 
(Atwood ct al. 2009) sensitive to 7 rays from 20 MeV 
to at least 300 GeV with on-axis effective area > 1 GeV 
of ~ 8000 cm^, exceeding that of EGRET by a factor of 
about five. It has a large field-of-view of 2.4 sr and in 
survey mode observes the entire sky every 3 hr. We re- 
port here on observations using data collected during the 
initial 35 days in on-orbit verification that included sky- 
survey tuning and pointed-mode tuning on Vela, from 
June 30 - August 3, 2008, as well as the initial 15 weeks 
of sky survey, from August 3 - November 16, 2008. The 
Diffuse class events (LAT event class having the tight- 
est background rejection) from these periods total 6014 
photons with energy > 100 MeV, within a radius of 1?5 
surrounding PSR J1028-5819. We excluded periods when 
the pulsar was viewed at zenith angle > 105° to the de- 
tector axis where the Earth's albedo photons gave exces- 
sive background contamination. 

Radio observations of PSR J1028— 5819 were carried 
out at the Parkes 64-m radio telescope at frequencies 
near 1.4 and 3.1 GHz with typical durations of ^5 min. 
Timing observations commenced on 2008 April 7, shortly 
after the discovery of the radio pulsar. Since then 24 in- 
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dependent timing measurements have been made with 
typical uncertainties in the times-of-arrival of ^30 /is or 
better. The fit to the timing points was carried out us- 
ing TEMP02 (Hobbs et al. 2006); the position of the 
pulsar was held fixed at the value given in Keith et al. 
(2008) and the rotation frequency and frequency deriva- 
tive were fit. The dominant contribution to the resulting 
residual of 270 /ks is timing noise intrinsic to the pulsar. 
Note that the timing solution used to derive the narrow 
radio profiles contains an extra fit parameter that is not 
included in the Fermi software tools used for the 7-ray 
profile. Including this additional parameter leads to at 
most a phase shift of 0.01 over the dataset used for the 
7-ray pulsation analysis. 

3. RESULTS 

For the detection of pulsations, Diffuse class events 
with energy > 100 MeV and within a radius of 1?5 of 
the radio position were corrected to the Solar-System 
barycenter using the JPL DE405 Solar System ephemeris 
and folded with the radio period using the Parkes 
ephemeris. The Fermi LAT timing is derived from a GPS 
clock on the satellite and photons are timestamped to an 
accuracy better than 300 ns. The LAT software tools 
for pulsars have been shown to be accurate to a few /US 
for isolated pulsars (Smith et al 2008). We detect 7-ray 
pulsations at the radio period with chance probability 
2 X 10^^^ using a Z„-test with 2 harmonics (DcJager et 
al. 1989). Within 2 weeks from the onset of data collec- 
tion, a Zn significance of 3.5 ct was found for pulsations, 
and this was improved to a better than 10 cr pulsed sig- 
nal with less than 4 months of LAT data. The 7-ray 
pulsations at the same period were independently found 
by a blind search in P and P using a time-differencing 
technique (Atwood et al. 2006, Ziegler et al. 2008). 

The LAT has an angular resolution with a dependence 
on reconstructed event energy E of E~'^-'^^, with a 68% 
containment radius of ~ 0? 5 at 1 GeV for near on axis 
civents which convert in one of the first 12 layers of the 
tracker (the thin section) and that increases with inci- 
dence angle as detailed in Atwood et al (2009). Events 
converting in one of the last four layers of the tracker 
(the thick section) have a 68% containment radius which 
is 2 times that of the thin section. In order to ex- 
plore the energy dependence of the light curve, the event 
selection was refined to be 0c{E/lOOMeV)~'^-'^^ , where 
0c = 3° for thin events and 9c ~ 4?1 for thick events 
are the containment radii at 100 MeV chosen to maxi- 
mize the detection significance with this energy depen- 
dent cut. The histogram of folded counts at energies 0.1 
-13 GeV is shown in Figure 1. The 7-ray light curve 
shows two strong peaks, PI at phase 0.200 ± 0.003 and 
P2 at phase 0.661 ±0.002, where phase is defined by the 
dedispersed radio pulse. The phase separation of PI and 
P2 is 0.460 ± 0.004. The peaks are fairly narrow, with 
Lorentzian FWHM for PI of 0.040 ± 0.011 and for P2 of 
0.035 ± 0.007. This light curve is very similar to that of 
Vela, as seen with EGRET (Kanbach et al. 1994) and 
now with Fermi (Abdo et al. 2009a). Figure 2 shows 
light curves in four different energy ranges, 100 — 300 
MeV, 300 MeV - 1 GeV, 1-3 GeV and > 3 GeV that do 
not exhibit significant evolution in shape. In particular, 
we have measured the widths of the peaks as a function 
of energy and for PI there is a deviation at only the 2.19 



a level and for P2 at the 1.63 a level. There is also no 
evidence for significant evolution of the P1/P2 ratio. A 
fit of the P1/P2 values to a constant function of energy 
gives a per degree of freedom of 0.54, which is consis- 
tent with no variation. This constrasts to the decrease 
in the P1/P2 ratio with energy seen in the Crab, Vela, 
Geminga and B1951+32 pulsars by EGRET (Thompson 
2004) , and is now confirmed in the Vela pulsar by Fermi 
(Abdo et al. 2009a). We have measured the pulsed sig- 
nificance of PSR J1028-I-5819 as a function of increasing 
energy, successively raising the low energy threshold. For 
all events above 4 GeV (with the energy dependent cut) 
the pulsed detection is 3.5 a, but above 5 GeV it is only 
1.8 cr, indicating that the maximum energy of pulsations 
is around 4 GeV. 

The LAT point source OFGL J1028.6-5817 from 
the Fermi LAT bright source list (Abdo et al. 
2009b) corresponding to PSR J1028-5819 is located 
at (RA,DEC)=(157.166,-58.292) with a 95% confidence 
level radius of 0.079 degrees. There are two other LAT 
point sources nearby, OFGL J1024.0-5754 and OFGL 
J1018.2-5858, 0?73 and 1?52 away respectively. The 
COS-B source 2CG 284-00 (Swanenburg et al. 1981) 
was apparently made up of contributions from all three 
LAT sources while the EGRET source 3EG J1027-5817 
has now been resolved by the Fermi LAT into contri- 
butions from the two sources. OFGL .11028.6-5817 and 
OFGL J1024.0-5754. The LAT source associated with 
PSR J1028-5819 is therefore somewhat confused with 
the nearby point sources and contains significant photon 
flux from OFGL J1024.0-5754. This adds to the unpulsed 
background and must be taken into account in computing 
the phase-averaged flux from the pulsar. A Swift XRT 
(Burrows et al. 2005) observation on November 23-24, 
2008 (9.6 ks cleaned exposure) yielded a 4.1cr detection 
of an X-ray source that we tentatively associate with the 
pulsar. Although the X-ray source is positionally coinci- 
dent with the radio pulsar, it is poorly localized (~20" 
total extent) and is faint, with a net count rate (0.3-10 
keV) of 2.3 X 10~^ cts/sec which corresponds to an ab- 
sorbed (unabsorbed) flux of 1.5 (3.4) x 10~^^ erg cm~^ 
s~^, assuming a power-law with r=2 and A^h = 1-59 x 
10^2 cm-2. 

To obtain the phase-averaged flux of the pulsar, we 
have performed a maximum likelihood spectral analysis 
using the LAT tool gtlike^'^ on counts within a radius 
of 15° from the LAT source position. The pulsar spec- 
trum was fit with a power law with an exponential cut- 
off, giving an index 1.22 ± 0.2 ± 0.12 and cutoff energy 
of 2.5 ± 0.6 ± 0.5 GeV (the first errors are statistical 
and the second are systematic). With the current statis- 
tics in this crowded region, a differentiation between a 
simple and a hyper-exponential cutoff is not yet possi- 
ble and, therefore, was not attempted. From this fit, 
we obtain an integral photon fiux at 0.1 — 30 GeV of 
1.62 ± 0.27 ± 0.32 x lO-'^phcm-^ s"^ and integral en- 
ergy fltix of 1.78 ±0.15 ±0.35 x lO'^^ergcm-^ s'^. This 
flux is a quarter of the EGRET source flux, which ap- 
parently included contributions from both PSR J1028- 
5819 and its neighboring source OFGL J1024.0-5754. 
In addition to the pulsar, the diffuse 7-ray emission 
from the Milky Way was modeled, while the extra- 

http: / / fermi.gsfc.nasa.gov/ssc/data/ analysis/SAE.overview.html 
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galactic diffuse emission plus instrumental residual back- 
ground and other LAT point sources within the region 
of interest were fit with power law spectra. The LAT 
source OFGL J1024. 0-5754 which likely contributed to 
the EGRET source was fit with an integrated flux from 
0.1 - 30 GeV of 3.0 x lO-'^phcm'^ s'^ The flux from 
this source plus that from PSR J1028-5819 add up to 
^ 4.8 X 10~^phcm~^ s~^ which is consistent, within 
statistical and systematic uncertainties, with the flux 
6.6 ± 0.7 X IQ-^phcm-^s-i of 3EG J1027-5817. The 
3EG source flux was derived by assuming a power law 
with index 2 (Hartman et al. 1999) which would have 
preferentially weighted the spectrum with a larger PSF, 
adding more counts from the background at low energy to 
produce a higher flux. The LAT source OFGL J1018.2- 
5858, which likely contributed to the flux of the 2CG 
source, yielded an integrated flux at 0.1 — 30 GeV of 
5.41 X 10~^ phcm~^ s~^ from a power-law fit. Adding 
the flux from PSR J1028-5819, OFGL J1024.0-5754 and 
OFGL J1018. 2-5858 then could account for the 2CG flux 
of ~ 2.7 X lO-'^phcm-^s-i. 

4. DISCUSSION 

The detection of pulsed 7 rays from the recently discov- 
ered pulsar PSR J1028-5819 within the first few weeks of 
the Fermi mission confirms the promise that the Fermi 
LAT will be an important instrument for pulsar stud- 
ies. The 7-ray pulses of PSR J1028-5819 cover a wide 
phase range and neither of the peaks is aligned with 
the very narrow radio pulse. This strongly suggests that 
the 7-ray beam covers a large fraction of solid angle of 
the sky and favors its interpretation in outer magneto- 
sphere models such as the outer gap (OG) (Cheng, Ho & 
Ruderman 1986, Romani & Yadigaroglu 1995) or slot 
gap (SG) (Mushmov & Harding 2004). Furthermore, 
the maximum observed energy of pulsations, tmax — 4 
GeV must lie below any 7-B pair production turnover 
threshold, thereby providing a lower bound to the alti- 
tude of emission, even though we have not been able to 
rule out the hyper-exponential spectral cutoff expected 
for pair production attenuation. Using a standard po- 
lar cap model estimate for the minimum emission height 
of r > (e™a^Bi2/1.76GeV)2/7p-i/7^^ (inverting Eq. 
[1] of Baring 2004), for a surface polar field strength 
of B12 = _B/10^^ G and neutron star radius i?*, the 
PSR J1028-5819 spin-down parameters (P = 0.0914s, 
B12 = 1.015) yield r > l.Si?,. This bound precludes 
emission very near the stellar surface. 

In order to estimate the 7-ray efficiency of a pul- 
sar, one needs to know the total luminosity radiated, 
L*°* = A-nf{aXE)^obsd^, where / is a correction factor 
that contains information about the beaming geometry, 
$of)s is the observed phase-averaged energy flux, and d 
is the distance to the source. The factor f{a,C,E) is a 
function of the pulsar magnetic inclination angle a and 
the viewing angle to the rotation axis, Ce, is very model 
sensitive and can be computed for any particular model 
geometry by, 

jF{aX,4')dn 



(1) 



(Walters et al. 2008). Here, F{aX,4') is the radiated 
fiux from a pulsar as a function of inclination angle, view- 
ing angle and rotation phase, 0. The factor /(a, (^b) is 



very important because it quantifies the amount of the 
emission we may be missing with our limited sweep over 
the pulsar beam. For polar cap models where the emis- 
sion originates within several stellar radii of the neutron 
star surface, the effective emission solid angle is small 
and thus the factor / <C 1. For outer magnetosphere 
models, the traditional solid angle measure is not appro- 
priate since the emission is radiated over a large fraction 
of An. The beaming correction factor / must therefore 
be computed numerically using Eqn ((T|) and it is found 

that / ^ 1 for both OG and SG models (Walters et 
al. 2008). For PSR J1028-5819, we obtain a total lumi- 
nosity of = 1.1 X lO'^^/ergs"^ from the observed 
energy flux and source distance, d ~ 2.3 kpc. The 7- 
ray efficiency is thus rjj = Ll!j'^/Esd ~ 0.13 ///45, where 
/45 — //lO'^^gcm^ is the neutron star moment of iner- 
tia. Since the distance could have as much as a 40% error 
from fluctuation in the free electron density (Brisken et 
al. 2002), there is an uncertainty of about a factor of 2 
in the derived luminosity and the efficiency. If the outer 
magnetosphere interpretation of the pulsed 7 rays from 
PSR J1028-5819 and other young pulsars is correct, then 
these pulsars have efficiencies that are larger by about 
an order of magnitude than would be deduced using the 
previously standard 1 sr solid angle. 

The fuU geometry of PSR J1028-5819 is not easily de- 
termined from the radio data alone because the narrow- 
ness of the pulse make it difficult to derive good solutions 
from polarization position angle variation using the ro- 
tating vector model. But the very narrow pulse does 
argue for a large a (in agreement with the constraints 
below). In the outer-magnetosphere geometry now fa- 
vored by the 7-ray emission, one can derive constraints 
on a and C, from the measured separation of the 7-ray 
peaks and 7-ray efficiency alone. Using the 7-ray light 
curve 'Atlas' of Walters et al. (2008), we estimate an al- 
lowed range of a ~ 70° - 90°, C ~ 75° - 80° and / ~ 1.1 
for the OG model and a - 65° - 80°, C ~ 60° - 80° and 
/ - 0.9 - 1.0 for the two-pole caustic (TPC) or SG mod- 
els. Both of these models therefore have a good range for 
viable solutions. However, the above estimates assume 
the 7-ray efficiency relation ry^ ~ (lO'^^ ergs~^/ii^id)^^^ 
which gives rj^ = 0.03, a factor of 4 smaller than that de- 
rived for this pulsar from its measured luminosity. Using 
our derived rj^ = 0.13, the range allowed for a shrinks 
to ~ 80° — 90° for the OG model, but remains about 
the same for the TPC model. Detection of a pulsar wind 
nebula (which can give an estimate of Cb) or radio polar- 
ization measurement (which might give an a estimate) 
would help further constrain the model. The phase lag 
of the first 7-ray peak relative to the radio pulse can be 
explained in either OG or SG (or generally TPC, Dyks 
& Rudak 2003) models, where the radio phase crossing 
associated with a magnetic pole occurs before the first 
high-energy emission caustic formed at high altitude. 

If many of the other young radio pulsars discovered in 
EGRET error boxes are similar to PSR J1028-5819 and 
if the 7-ray beams of young pulsars can be seen from 
a wide range of viewing angles, then Fermi searches for 
their pulsed 7 rays promise to be very fruitful. The idea 
of wide 7-ray beams recently gained strong support from 
the Fermi LAT discovery of a radio quiet pulsar in the 
young supernova remnant CTA 1 (Abdo et al. 2008) 
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whose pulsations have been detected only in 7 rays, a 
blind search result secured promptly during the Fermi 
commissioning phase. The discovery of 7-ray pulsations 
from a young radio pulsar coincident with an EGRET 
source, together with the discovery of 7-ray pulsations 
from CTA 1, confirms expectations from before the Fermi 
launch that many 7-ray pulsars remain to be discovered. 
Even more exciting are the prospects that the all-sky 
survey of the Fermi LAT will identify many new sources, 
and that searches in both the 7-ray and radio bands will 
uncover an entirely new population of pulsars. 
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Fig. 1. — Light curve of PSR J1028-5819 in the (0.1 - 13 GeV) band in 40 constant-width bins and shown over two pulse periods with 
the 1.4 GHz radio pulse profile plotted below. The horizontal dashed line shows the estimated background level from the off-pulse region 
at phase 0.8 and 1.0. The inset shows the radio pulse in the phase range 0.97 - 1.03, with the main peak at phase 1.0 and preceded by the 
smaller, secondary peak at phase ~ 0.996. 
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Fig. 2. — Light curves of PSR J1028-5819 in four different energy bands (labeled) in constant width bins of size 0.04 in phase. The 
horizontal dashed lines show the estimated background level from the off-pulse region at phase 0.8 and 1.0. 



